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Abstract 

We propose a new method to extract the light quark mass ratio m u /rrid using the T(4S) — > hvK \T\) bottomonia transitions. The decay 
amplitudes are dominated by the light quark mass differences, and the corrections from other effects are rather small, allowing for a precise 
extraction. We also discuss how to reduce the theoretical uncertainty with the help of future experiments. As a by-product, we show that the 
decay T(4S f ) —¥ hbrj is expected to be a nice channel for searching for the hb state. 



PACS numbers: 13.25.Gv, 14.65.Bt, 12.39.Hg 

Although fundamental parameters of the Standard Model, 
the masses of light quarks have not yet been well determined. 
This appears to be a consequence of quark confinement as 
well as the fact that the light quark masses are significantly 
lighter than the typical hadronic scale and as such their im- 
pact on most of the hadron masses or other properties is very 
small. 

As a consequence of the spontaneous chiral symmetry 
breaking, the low-energy region of the quantum chromody- 
namics (QCD) can be described by chiral perturbation theory 
(CHPT) |ll |2y. The most direct way to get information on 
the light quark mass ratios is to relate the quark masses to 
the masses of the lowest-lying pseudoscalar mesons, which 
are the Goldstone bosons of the spontaneously broken chi- 
ral symmetry of QCD. To leading order (LO) in the chiral 
expansion, this gives m u /md = 0.56 yfl. Electromagnetic 
(e.m.) effects have been taken into account using Dashen's 
theorem ]4[]. There might be, however, sizeable higher or- 
der corrections to this LO result, e.g. related to violations of 
Dashen's theorem, see |5|,|g]. The up-to-date knowledge of the 
light quark mass ratio from various sources including recent 
lattice calculations was summarized in Ref. |7J] to be 
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m d 



0.47 ±0.08 



(1) 



In a completely independent approach it was proposed to 
use the decays of ijj' into J/ipir° and J jtbr\, which break 
isospin and SU(3) symmetry, respectively 181I9J. It was as- 
sumed that these decays are dominated by the emission of soft 
gluons, and the gluons then hadronize into a pion or an eta. 
Using the QCD multipole expansion (QCDME), one obtains 
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where F n ^ and M„(rf\ are the decay constant and mass of 
the pion (eta), respectively, and q^M is me pi° n ( eta ) momen- 
tum in the ip' rest frame. These two decays were widely used 
in determining the quark mass ratio m u /md II 1 OJ - 4 1 311 . Using 
the most recent measurement of the decay widths from the 
CLEO Collaboration 0], one gets m u /rn d = 0.40 ± 0.01, 
which is much smaller than the one resulting from the meson 
masses. Using instead the measurement by the BES Collab- 
oration || 1 5-1 , the resulting value m u /md — 0.35 ± 0.02 is 
even smaller. In Ref. jlql , based on a non-relativistic effective 
field theory (NREFT) formalism, the striking discrepancy be- 
tween the values of m u /md extracted from the ij/ decays and 
from the meson masses was solved by showing that the decay 
amplitudes of the transitions ip' —> J/tpir (77) are not domi- 
nated by the multipole effect as assumed before. Rather, non- 
multipole effects via intermediate charmed meson loops are 
very important, enhanced by 1/v, V ~ 0.5 being the charmed 
meson velocity, compared with the multipole one. More pre- 
cisely, the large uncertainty related to the non-multipole con- 
tributions prevents one from an extraction of m u /md from 
these decays. 

In this Letter, we propose a new way to extract the light 
quark mass ratio using the transitions of the excited bottomo- 
nium T(45) into hhir and hi,?]. Similar to the transitions be- 
tween charmonium states, the e.m. contribution to the isospin 
breaking decay T(4S) — > /ifc7r° is negligibly small Il7l4l9ll . 
This provides the possibility of extracting the light quark 
mass ratio from these decays. It will be shown that the non- 
multipole effects from intermediate bottom meson loops are 
suppressed, and hence the decay amplitudes are proportional 
to the light quark mass differences. 

The spin-singlet P-wave bottomonium hb has not been ob- 
served yet, however, it is expected to agree in mass with the 
spin-averaged mass of the spin-triplet P-wave bottomonia \b j 
(see, e.g., Ref. Q), which is M hb = 9900 MeV. The T(45) 
with a mass of 10579.4 ± 1 .2 MeV andwidth 20.5 ± 2.5 MeV 
is the first bottomonium above the BB threshold, and it de- 
cays into BB with more than 96% branching fraction 12 ill . 
The mass difference between the T(45) and the hb is about 



680 MeV. Hence, both the transitions T{4S) -t h b ir and 
T(4S) — ?> h b T] are kinematically allowed. 

Let us consider the multipole decay mechanism with the 
light meson being directly emitted from the bottomonium 
first, which is described by a tree-level diagram based on 
hadronic degrees of freedom. Because the decays are in an 
S-wave and break isospin or SU(3) symmetry, the LO ampli- 
tude must scale as the quark mass difference 
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(3) 



with 5 = m,d — m u for the transition T(45) — > h b ir and 
5 = m s — rh, with rh = (m u + m<j)/2, for the transition 
T(45) -> h bV . 

Corrections to the tree-level result arise due to intermedi- 
ate heavy meson loops and higher order terms in the chi- 
ral expansion. The loops can be studied in the framework 
of the NREFT because the velocity of the heavy meson in 
the loops is small. The value of the bottom meson veloc- 
ity for the transitions considered here may be estimated as 
v ~ y/[2M 6 - (M T(4S) + M hb )/2\/M 6 ~ 0.3 with M 6 
the averaged bottom meson mass. This estimate is consis- 
tent with determinations of the bottom quark velocity in bot- 
tomonium systems based on non-relativistic QCD (see, e.g. 
Ref. IJ22ID . For a transition between a P-wave and an S-wave 
heavy quarkonium with the emission of a pion or an eta, it has 
been shown that the contribution to the decay amplitude from 
the intermediate heavy meson loops scales as 12311 
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H 



where q is the three-momentum of the light meson in the rest 
frame of the decaying heavy quarkonium, Mg is the mass of 
the intermediate heavy meson, and the meson mass difference 
A encodes the violation of the isospin symmetry for the pio- 
nic transition or SU(3) symmetry for the eta transition. Eq. (HJi 
arises because the non-relativistic loop integral measure con- 
tains three powers of momentum, and scales as v 3 . After per- 
forming the contour integration of the energy, two propagators 
are left, and each of them scales as 1/v 2 . The P-wave cou- 
pling of the light meson to the heavy meson gives a factor of q. 
The coupling of the heavy mesons to the P- and S-wave heavy 
quarkonia are in S- and P-wave, respectively. The P-wave ver- 
tex provides a momentum in the loop integral, and it must be 
contracted with the external momentum of the light meson q. 
So the three vertices together provide a factor of q 2 . Since 
we are considering isospin or SU(3) symmetry breaking tran- 
sitions, the decay amplitude from the loops is non-vanishing 
because the heavy mesons within the same isospin or SU(3) 
multiplet have different masses. One may pull out the meson 
mass difference explicitly to represent the symmetry break- 
ing. Because it is an energy scale and should be counted as 
v 2 , one needs to divide it by v 2 for balance. Putting all pieces 
together, one gets [v 3 /(v 2 ) 2 ][q 2 /M H -][A/v 2 ], where 1/Mj|- 
is introduced to match dimensions, and Eq. (@]i follows. This 
kind of non-relativistic power counting has already been con- 
firmed by explicit calculations of the loops jig, ll9L 12311 . 



To determine the relative size of the loop amplitude com- 
pared to the tree-level one, in order to find out whether 
the tree-level contribution is dominant, one should compare 
the meson mass difference A and the quark mass differ- 
ence 8, and estimate the value of the dimensionless pre-factor 

2 /(v 3 Mff). The momenta of the pion and eta in the final 



9 



states of T(4S) -> h b ir°(v) are 645 MeV ad 389 MeV, re- 
spectively. Taking v w 0.3 for the velocity, the dimensionless 
factor q 2 /(v 3 M B ) is about 0.6 for the pionic transition and 
0.2 for the eta transition. One cannot naively assign the me- 
son mass differences as the same order as the quark ones. In 
fact, due to destructive interference between the quark mass 
difference and the e.m. contribution 12411 . the isospin mass 
splitting of the bottom mesons B° and B + is rather small, 
M B o-M B + = 0.33 ±0.06 MeV J2H. It is one order of mag- 
nitude smaller than md—rn u . Together with the dimensionless 
factor, which is about 0.6, the bottom meson loops contribute 
to the decay T(4S) — > hbir for no more than a few percent, 
and hence are negligible. The situation for the eta transition is 
somewhat different because M Bs - M B = 87.0 ± 0.6 MeV, 
where M B = (M B o + Mg+)/2, and it is of similar size 
as m s — rh. This means that the loop contributions to the 
T(45) — ^ h b r\ as compared to the tree-level decay amplitude 
are also suppressed, but they might give a non-vanishing cor- 
rection of about 20%. 

An intriguing implication of the suppression of the bottom 
meson loops in these transitions is that the decay amplitudes 
are dominated by the quark mass differences, and hence it is 
possible to extract the light quark mass ratio from the ratio of 
the branching fractions of the transitions T(4S) — > h b 7r°(r)) 
with good accuracy. It has been demonstrated that the LO 
results of chiral Lagrangians for the heavy quarkonia transi- 
tions can reproduce the LO results of the QCDME 12511 . In 
the QCDME, the transitions between two heavy quarkonia 
occur through radiating soft gluons, and the soft gluons then 



2811 . In the case of transi- 



hadronize into light mesons (26 
tions with the emission of a pion or an eta, the gluon operator 
is GG = a s G^ v G A ^ v H, where a s is the strong coupling 
constant, G A V is the gluon field strength tensor and its dual is 
QAiiu = e wc£ a /2. For the transitions T(4S) -> h b ir (r)), 
we have 



r(T(4S) -> h b ^) 
r(T(45) -► h b rj) 



GG 



(5) 



where q^uj) is the momentum of the pion (eta) in the rest 
frame of the T(45), and the ratio of the gluon matrix elements 
is defined as 



' GG 



(0|GG|7T°) 
(0\GG\ v ) 



(6) 



Combining CHPT with the U(1)a anomaly, the next-to- 
leading order (NLO) expressions for the matrix elements 
(0|GG|7r°(?7)) were worked out in Ref. [1311 m terms of sev- 
eral low-energy constants (LECs) of the 0(p 4 ) Lagrangian. 
Moreover, there exists an intriguing relation between the ratio 



of the matrix elements and a combination of the light quark 
masses II 1 1 L 1 1 3fl 
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= 10.59(1 + 132.1 L 



14J 



GG> 



(7) 



where <5gmo = —0.06 denotes the 0(p A ) deviation from 
the Gell-Mann-Okubo relation among the Goldstone bosons. 
Higher order terms in the coupling of the flavor-singlet field 
that encodes the information of the anomalously broken 
U(1)a anomaly are parameterized by the 0(p 4 ) LEC L14. 
Therefore, once one has knowledge of the value of the LEC 
L14, one is able to extract the quark mass ratio from the ra- 
tio of branching fractions of the decays T(4S) — > h b ir° and 
T(4S f ) — > hbT), which can be measured in the future. 

There are two main theoretical uncertainties for extracting 
the value of row ■ The first one is due to lack of knowledge 
of the LEC L14. One may use resonance saturation to esti- 
mate its value, and it is expected to be in the region llll 1 1 3fl 
L u = (2.3 ± 1.1) x 10~ 3 . From Eq. 0, it gives 11% un- 
certainty in ?*dw- The other one is from neglecting the inter- 
mediate bottom meson loops of the transition T(4S) —> h b rj. 
As already discussed, it gives an uncertainty of 20% in the 
amplitude, and hence 40% in the decay width. Propagating to 
the extracted quark mass ratio, the uncertainty is again 20%. 
Adding them quadratically, the theoretical uncertainty for ex- 
tracting tdw is 23% which is comparable to that of Eq. (H). 

The uncertainty could be reduced once further information 
on the size of the loops is available. This kind of information 
could be provided by high statistics measurements in the fol- 
lowing way: The decay width for T(45) — i h b rj considering 
only loops can be worked out using the NREFT 



T(T(4S) 



M)'° op = 0.16^ fc keV, 



(8) 



where the only unknown parameter g\ b denotes the coupling 
of the IP bottomonium states to the bottom mesons, given 
in GeV -1//2 . Although the IP states are below any open 
bottom threshold, one may extract g lb from the loop domi- 
nated transitions involving the IP states. Because the isospin 
mass splitting of the B mesons is rather small, one should 
consider the loop dominated transitions with the emission of 
an r\. These are the transitions from excited P-wave bot- 
tomonia to the IP states, enhanced by a factor 1/v 3 , using 
a similar power counting technique presented in Ref. 11911 . 
The best choice are the r\ transitions from the 4P states to 
the IP states. Based on the quark model calculation of the 
bottomonium spectrum 02911 . the 4P states have sufficiently 
large masses to allow for decays into B^*'&-*> , Because the 
4P states can decay directly into B^B^*\ the corresponding 
coupling constant g'^ can be obtained by measuring their de- 
cay widths. Then, one can extract the value of gu, from any 



of the transitions x &0 (4P) ->• XbiV> Xfei( 4 -P) -► X60,i,2»7, 
Xm(4P) -¥ Xbi,2V and /i b (4P) -> h b r\. In Fig. [TJ we 
show the predictions from the NREFT for the following ra- 
tios, which depend on gu, only and are proportional to g\ h , as 
a function of the mass of the 4P bottomonium state: 



P01 



R\j= 



R2J 



r( Xfc0 (4P) -> xbiv) 

T{ Xw ^P)^B+B-y 
r( Xb i(4P) ->■ xbjrj) 

r( Xfc i(4P)^p+p*-)' 

r(x b2 (4P) -> xbjy) 
r(xb 2 (4P)^p*+p*-)' 



J = 0,1,2, 



J = 1,2. 



(9) 



The result for T(h b (4P) -> h b r))/T(h b (4P) -)■ B*+B*-) 
is very similar to and slightly larger than Pin- The cusps in 
Fig.[TJ(a) and (b) represent the opening of the B S B S and B S B* 
thresholds, respectively. For definiteness, we have used g± b = 
1 GeV" 1 / 2 . The dependence on g'H is canceled in the ratios. 
If any of these ratios were to be measured, one will be able to 
extract the value of g\ h easily. The so extracted coupling g\ h 
bears about 30% uncertainty due to the loops of higher order. 
The 4P bottomonia should decay dominantly into B^B^*\ 
For an order-of-magnitude estimate, one may assume g lb to 
have a similar value as its charm analogue estimated using 
vector meson dominance, <? lc = —4.2 GeV -1 / 2 Q30Q . Then 
from Fig. [TJ one expects the eta transitions have branching 
fractions of the order of a few per cent. Hence, there should 
be a good opportunity in extracting g\ b at the Large Hadron 
Collider beauty experiment (LHCb). After having measured 
the partial decay width of T(4S) — > h b i], one may compare 
the measured value with the one obtained considering only the 
bottom meson loops, given in Eq. ((H), using gi b determined in 
the way outlined above as input. Depending on whether the 
interference between the tree-level and the loop amplitudes is 
constructive or destructive, one can get two solutions of the 
width considering only the multipole (tree-level) effect. Then 
one may insert the resulting T(T(45) —¥ h b r]) trc ° in Eq. (0, 
reducing the uncertainty from the loop contribution. 

Using the same naturalness arguments for the coupling con- 
stant in the LO tree-level Lagrangian as that in Ref. 12311 . 
we can estimate the branching fractions of the transitions 
T(45) — > h b TT°(r]). The branching fraction for the pionic 
transition is of order 10 -6 , and the one for the eta transition 
is of order 10 -3 . With such a large branching fraction, the 
latter one even provides a nice option for searching for the h b . 
LHCb is expected to have enough events of the T(4S I ) to do 
the measurements. 

In summary, we have proposed a new method for extract- 
ing m u /m,d- We demonstrated that the transitions T(4S') — > 
h b Tr°(r]) can be used to determine the value of row with an 
acceptable theoretical uncertainty, which is about 23%. Us- 
ing information of rh/m s from other sources, one is then 
able to extract m u /md- The transitions T(45) — > h b ir° and 
T(45) — > h b rj are expected to have branching fractions of or- 
der of 10~ 6 and 10 -3 , respectively. Therefore, they can be 
measured at LHCb based on a large number of T(4S f ) events. 
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FIG. 1 : Predicted ratios defined in Eq. ((9} in the NREFT. The value of gn, has been set to 1 GeV 1 ' 2 . 



The uncertainty can be reduced to obtain a more accurate ex- 
traction of the quark mass ratio by measuring the partial decay 
widths of the 4P bottomonium to the IP bottomonium with 
the emission of an eta. These transitions with branching frac- 
tions of order of a few per cent can also be measured at LHCb. 
As a by-product, the decay T(4S) — > h\,rj is expected to be a 
nice channel for searching for the hb state. 
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